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1 Introduction 



Although the main subject of this article is the connection between integrable dis- 
crete systems and geometry, we feel obliged to begin with the differential part of the 
relation. 



1.1 Classical differential geometry and integrable systems 

The oldest (1840) integrable nonlinear partial differential equation recorded in liter- 
ature is the Lame system 



d 2 Hi 1 dHjdH, 1 8H k dHi 



dujduk Hj duk duj Hk duj duk 



0, k distinct (1.1) 



d ( 1 dHj\ d ( 1 dH k \ 1 dHjdHk . . 

+ 772^^ = °> ( L2 ) 



duk \Hk duk ) duj \Hj duj ) Hf dui du 

describing orthogonal coordinates in the three dimensional Euclidean space E 3 (in- 
dices k range from 1 to 3). Already in 1869 it was found by Ribaucour that the 
nonlinear Lame system possesses a discrete symmetry enabling to construct, in a 
linear way, new solutions of the system from the old ones. He gave also a geometric 
interpretation of this symmetry in terms of certain spheres tangent to the coordinate 
surfaces of the triply orthogonal system. In 1918 Bianchi showed that the result of 
superposition of the Ribaucour transformations is, in a certain sense, independ of the 
order of their composition. 

Such properties of a nonlinear equation are hallmarks of its integrability, and in- 
deed, the Lame system was solved using soliton techniques in 1997-98. The above 
example illustrates the close connection between the modern theory of integrable 
partial differential equations and the differential geometry of the turn of the XlXth 
and XXth centuries. A remarkable property of certain parametrized submanifolds 
(and then of the corresponding equations) studied that time is that they allow for 
transformations which exhibit the so called Bianchi permutability property. Such 
transformations called, depending on the context, the Darboux, Calapso, Christoffel, 
Bianchi, Backlund, Laplace, Koenigs, Moutard, Combescure, Levy, Goursat, Ribau- 
cour or the fundamental transformation of Jonas, can be geometrically described in 
terms of certain families of lines called line congruences. 

In the connection between integrable systems and differential geometry, a dis- 
tinguished role is played by the multidimensional conjugate nets, described by the 
Darboux system, which is just the first part (jl.lj) of the Lame system with indices 
ranging form 1 to iV > 3. On the level of integrable systems, this dominant role 
has the following explanation: the Darboux system, together with equations describ- 
ing iso-conjugate deformations of the net, forms the multicomponent Kadomtsev- 
Petviashvilii (KP) hierarchy, which is viewed as a master system of equations in soli- 
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ton theory. In fact, in appropriate variables, the whole multicomponent KP hierarchy 
can be rewritten as an infinite system of the Darboux equations. 

1.2 Transition to the discrete domain 

The recent progress in studying discrete integrable systems showed that, in many 
respects, they should be considered as more fundamental than their differential coun- 
terparts. Consequently, the natural problem of extending the geometric interpretation 
of integrable partial differential equations to the discrete domain arose, leading not 
only to the transition to the discrete domain of many results on the connection be- 
tween the differential geometry and integrable systems, but also - and this seems to 
be even more important - to the description of integrability in a very elementary and 
purely geometric way. 

On the level of integrable equations, the transition 'from differential to discrete' 
often makes formulas more complicated and longer. On the contrary, on the geo- 
metric level, in such a transition the properties of discrete submanifolds, relevant to 
their integrability, become simpler and more transparent. Indeed, the mathematics 
necessary to understand the basic ideas of the integrable discrete geometry does not 
exceed the 'ruler and compass constructions', and many proofs can be performed 
using elementary incidence geometry. 

We will concentrate our attention on the multidimensional lattice made from pla- 
nar quadrilaterals, which is the discrete analogue of a conjugate net. Together with 
the discussion of its properties, which are the core of the geometric integrability, we 
briefly present the analytic methods of construction of these lattices and we also de- 
scribe some basic multidimensional integrable reductions of them. Then we discuss 
integrable discrete surfaces; some of them have been found in the early period of 
the 'case by case' studies. We shall however try to present them, from a unifying 
perspective, as reductions of the quadrilateral lattice. 

2 Multidimensional integrable lattices 
2.1 The quadrilateral lattice 

An iV dimensional lattice x : Z N — > M M is a lattice made from planar quadri- 
laterals, or a quadrilateral lattice (QL) in short, if its elementary quadrilaterals 
{x, TiX, TjX, TiTjx} are planar; i.e., iff the following system of discrete Laplace equa- 
tions is satisfied 

AiA jX = (TiAij) AiX + (TjAjJAjX, i^j, i,j — l,...,N, (2.1) 

where Aij : Z N — > R are functions of the discrete variable; here 7$ is the translation 
operator in the ith direction and Aj = Tj — 1 is the corresponding difference opera- 
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Figure 1: The geometric integrability scheme 



tor. For simplicity we work here in the affine setting neglecting projective geometric 
aspects of the theory. 

2.1.1 The geometric integrability scheme 

In the case N = 2 the definition (j2.1)l allows one to uniquely construct, given two 
discrete curves intersecting in a common vertex and two functions A12, A21 '■ I* 2 — > 
R, a quarilateral surface. For N > 2 the planarity constraints ()2.1|) are instead 
compatible if and only if the geometric data satisfy the nonlinear system 



A k Aij + (TkAi^Aik = (TjAj^Aij + (T k A kj )A ik , i,j, k distinct. (2.2) 



This constraint has very simple interpretation: in building the elementary cube (see 
Figure HJ, the seven points x, TiX, TjX, T k x, TiTjX, T{T k x and TjT k x (i,j,k are 
distinct) determine the eighth point T{FjT k x as the unique intersection of three planes 
in the three dimensional space. 

The connection of this elementary geometric point of view with the classical theory 
of integrable systems is transparent: the planarity constraint corresponds to the set of 
linear spectral problems (|2.1|) and the resulting QL is characterized by the nonlinear 
equations ()2.2|) . arising as the compatibility conditions for such spectral problems. 
Since the QL equations ()2.2|) are a master system in the theory of integrable equations, 
planarity can be viewed as the elementary geometric root of integrability. The idea 
that integrability be associated with the consistency of a geometric (and/or algebraic) 
property when increasing the dimensionality of the system is recurrent in the theory 
of integrable systems. 
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2.1.2 Other forms of the Darboux system 

The % <-> j symmetry of the RHS of equations (|2.2|) implies the existence of the 
potentials Hi : Z N — > R (the Lame coefficients) such that 

^ = -5^ i±h (2-3) 
Hi 

and then equations ()2.2|) take the form 

A fc A^ - A 3 H t - [Tk^^-J A k H t = 0, k distinct, (2.4) 

being the discrete version of the first part (jl.lj) of the Lame system. 

The Lame coefficients allow to define the suitably normalized tangent vectors 
X,i : Z N — > R A/ by equations 

A lX = {TiHi)X i} (2.5) 
and the functions Qij : 7L N — > R, z 7^ j, (the rotation coefficients) by equations 

AiHj = (TiHjQn, i / ./. (2.6) 

Then equations (|2.1jl and (|2.4jl can be rewritten in the first order form 

A j X i = (T j Q ij )X j , i / ./. (2.7) 
AkQij = {TkQik)Qkj, i,j,k distinct. (2.8) 

The discrete Darboux system ()2.8j) implies the existence of other potentials pi 
defined by the compatible equations 



T jP i 



1 o;q j ,)(t j q,:). 1 / j. (2.9) 



Pi 

The i <-> j symmetry of the RHS of equations ()2.9|) implies existence of yet another 
potential r : Z N -> R, 

A = — , (2.10) 
T 

which is called the r-function of the quadrilateral lattice. In terms of the r function, 
and of the functions 

Tij = TQij, i^j, (2.11) 

whose geometric interpretation will be given in Section 12.3.31 the discrete Darboux 
equations take the following Hirota-type form 

(TiTjT)r = (TiT)TjT - (TiTji)TjTij, i ^ j, (2.12) 
(T k Tij)T = (T k T)Tij + (T k r ik )r kj , i,j,k distinct. (2.13) 
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2.2 Analytic methods 

We will show how one can construct large classes of solutions of the discrete Dar- 
boux equations and the corresponding quadrilateral lattices using two basic analyti- 
cal methods of the soliton theory: the d dressing method and the algebro-geometric 
techniques. 

2.2.1 The d dressing method 

Consider the non-local 9-problem 

B X (z) + (R X )(z) = 3v(z), lim ( X (z) - v(z)) = 0, (2.14) 

\z\^oo 

where 3 = d/dz, R is the integral operator 
(R X )(z) = [ R(z,z') X (z')dz' Adz', 



and u(z) is a given rational function of z. 

Let Qf e C, i = 1, . . . , iV be pairs of distinct points of the complex plane, which 
define the dependence of the kernel R on the discrete variable n G % N 



i=1 \ Tl / j = l 



We consider only kernels Rq(z,z') such that the non-local <9-problem is uniquely 
solvable. If x( z 'i n ) is the unique solution with the canonical normalization v — 1, 
then the function 



^(z-,n) = x(z-,n)f[( Z v -^] 

i=l ^ ^ l ' 



satisfies the system of the Laplace equations (j2.1j) with the Lame coefficients given 
by 

Hi(n)= lim ((^ppi \{z-n) 

-^Qi \ \ z -Qi J 

By construction, the system of such Laplace equations is compatible, therefore the 
Lame coefficients satisfy equations (|2.4j) . To various n-independent measures dfj, a on 
C there correspond coordinates 

x a (n) = / ip(z;n)djj a (z), 
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of a quadrilateral lattice x, having Hi(n) as the Lame coefficients. To have real 
lattices, the kernel R , the points Qf and the measures dfx a should satisfy certain 
additional conditions. 

One can find a similar interpretation of the normalized tangent vectors Xi and 
of the rotation coefficients Qij. If Xi( z 'i n ) are the unique solutions of the non-local 
(9-problem 1)2.14)1 with the normalizations 



then the functions ipi(z; n), defined by 

h— 1 ^ *fo / 



k=l 

satisfy the direct analogue of the linear problem (J2.7j) 

Ajipi(z;n) = (TjQijfaMjfan), i ^ j, (2.15) 

where 

Qij(n) = lim (( ^ I ipi(z;n) 



Again, by construction, equations ()2.15|) are compatible and the functions Q 



n 



satisfy the discrete Darboux equations (|2.8jl . The functions 

x i( n ) = / ipi(z;n)d/i a (z)i 



c 



are coordinates of the normalized tangent vectors Xi of the quadrilateral lattice x 
constructed above. 



2.2.2 The algebro-geometric techniques 

Given a compact Riemann surface 1Z of genus g, consider a non-special divisor 
D = J2a=i Pol- Choose N pairs of points Qf G 1Z and the normalization point 
Qoo. Given n 6 Z^, there exists a unique Baker- Akhiezer function ip(n), defined 
as a meromorphic function on 7Z, with the following analytical properties: (i) as a 
function of P G 7£\ U^ 1 Q^ t , ?/>(n) may have as singularities only simple poles in the 
points of the divisor D; (ii) in the points Qf function ip(n) has poles of the order 
±n«; (iii) in the point Qoo function ip(n) is normalized to 1. 
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When zf{P) is a local coordinate on 1Z centered at Qf, then condition (ii) implies 
that the function ip(n) in a neigbourhood of the point Qf is of the form 



i,{P;n) = {zt(P)f ni ££,±(n) • (2-16) 



,s=0 



The Baker- Akhiezer function, as a function of the discrete variable n G Z , satisfies 
the system of Laplace equations (jUT)) with the Lame coefficients i?j(n) = Co+( n )- 

Again, by construction, the Lame coefficients satisfy equations (|2.4j) . To various 
n-independent measures <i// a on 1Z there correspond coordinates 

x a {n) = [ V(P;n)rf// a (P), 



</7e 

of a quadrilateral lattice cc. 

We present the expression of the Baker-Akhiezer function and of the r-function 
of the quadrilateral lattice in terms of the Riemann theta functions. Let us choose on 
1Z the canonical basis of cycles {ai, . . . ,a g ,bi, . . . , b g } and the dual basis {a>i, . . . , oj g } 
of holomorphic differentials on 7Z, i.e., § a .^k — (^.Then the matrix B of 6-periods 

defined as Bjk = § b , tUfc,is symmetric and has positively defined imaginary part. De- 
note by ujpq the unique differential holomorphic in 1Z \ {P, Q} with poles of the first 
order in P, Q and residues, correspondingly, 1 and —1, which is normalized by con- 
ditions § a ujpq = O.The Riemann function 9(z; B), z G C 9 , is defined by its Fourier 
expansion 

9(z; B) = exp {ni(m, Bra) + 2ni(m, z)} , 

where (•, •) denotes the standard bilinear form in C 9 . Finally, the Abel map A is given 
by A(P) = y fp o u>i, . . . , fp o ujg^j , where P G 7Z, and the Riemann constants vector K 
is given by 

The explicit form of the vacuum Baker-Akhiezer function ip can be written down 
with the help of the theta functions as follows 



(a(P) + TL^k (A(Qk) - A(QD) + z) 

*p{n,P) = —± 

[A{Q oo) + Ef=i^ (A(Q- k ) - A(Qt)) + Z 



6 QA(goo) + Z) 
X B(A(P) + Z) 6XP 
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where Z = - A(Pj) - K. 



Denote by r^- and s k j the constants in the decomposition of the Abelian integrals 
near the point 



/ wg-gj P = Qj T5 kj log^(P) + r± + 0(*f(P)), 

JPo 

f " Qoo Qi P =^ "M+± l0g^(P) + ^ + 



Then the expression of the r-function of the quadrilateral lattice within the subclass 
of algebro-geometric solutions reads 

(N \ N N 

(A(Qi) - A(Ql)) + A(Q 00 ) + Z J J] A^ JI^' 
fc=i / fcj=i fc=i 

where 

Finally, we remark that the geometric integrability scheme and the algebro-geometric 
methods work also in the finite fields setting, giving solutions of the corresponding 
integrable cellular automata. 



2.3 The Darboux-type transformations 

We present the basic ideas and results of the theory of the Darboux type transforma- 
tions of the multidimensional quadrilateral lattice. 



2.3.1 Line congruences and the fundamental transformation 

To define the transformations we need to define first A^-dimensional line congruences 
(or, simply, congruences), which are familes of lines in R A/ labelled by points of Z, N 
with the property that any two neighbouring lines 1 and Til, i = 1, N, are coplanar 
and therefore (eventually in the projective extension P M of R A/ ) intersect. 

The quadrilateral lattice J-{x) is a fundamental transform of the quadrilateral 
lattice x if the lines connecting the corresponding points of the lattices form a con- 
gruence. The superposition of a number of fundamental transformations can be com- 
pactly formulated in the vectorial fundamental transformation. The data of the vec- 
torial fundamental transformation are: (i) the solution Y", : Z* N — > V, V being a linear 
space, of the linear system (|2.7|) : (ii) the solution Y* : Z N — > V*, V* being the dual 
of V, of the linear system ()2.6|) . These allow to construct the linear operator valued 
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Figure 2: The fundamental transformation as the binary transformation 
potential fl(Y, Y*) : Z N — > L(Y), defined by the following analogue of equation (|2.5jl 



A,n(Y, Y*) = Y t ® (TjY*), i = l,...,N; (2.17) 

similarly, one defines Q(X,Y*) : Z N -> L(V,M A/ ) and ft(Y", H) : -> V. The 
transforms of the lattice x and other related functions are given by 

T{x) = x - n(x, r*)rj(r, r^nor, 

JF(F 4 ) = Hi- Y*ft(Y, YT'tyY, H), i = l,...,N, 

r{Xi) = x, - n(x, r*)fi(y, r*)- 1 ^, « = 1, . . . , n, 

F(Qii) = Qij-Y*n(Y,Y*)- 1 Y h i,j = l,...,N, i^j, 
J={ Pl ) = Pi {l + (TiY^ftiY, Y*)Yi) , i = l,...,N, 
JF(r) = Tdettl(Y,Y*). 

Notice that, by the coplanarity of any two neighbouring lines of the congruence, also 
the quadrilaterals {x, TiX, !F(x), T{TiX)} are planar (see Figure EJ). Then the con- 
struction of the transformed lattice mimics the geometric integrability scheme. In 
consequence, any quadrilateral {x, ^(x), Ti{x), T\ (JF 2 (a;)) = JF 2 (^(a;))} is planar 
as well. Therefore, on the discrete level, there is no difference between the lattice co- 
ordinate directions and the fundamental transformation directions. The distinction 
becomes visible in the limit from the quadrilateral lattice to the conjugate net. There- 
fore the vectorial description of the superposition of the fundamental transformations 
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not only implies their permutability but also provides the explanation of the validity 
of the practical rule of "integrable discretization by Darboux transformations". 



2.3.2 The Levy and Combescure transformations 

It is easy to see that the family tj of lines passing through the points x and %x of 
a quadrilateral lattice forms a congruence, called the i-th tangent congruence of the 
lattice. When the congruence of the transformation is the z-th tangent congruence of 
the lattice x, then the corresponding reduction of the fundamental transformation is 
called the Levy transformation Li. 

It turns out that, for a generic congruence [, the lattice made from intersection 
points of the lines [ and T~ x \ is a quadrilateral lattice, called the i-th focal latice of 
the congruence. When the fundamental transform of the lattice x is the z-th focal 
lattice of the transformation congruence, then the corresponding reduction of the 
fundamental transformation is called the adjoint Levy transformation £*. 

Both Levy transformations use only a half of the fundamental transformation 
data, and the corresponding reduction formulas (in the scalar case) for the lattice 
points read as follows 



Notice that the composition of the Levy and the adjoint Levy transformations gives 
(see Figure 12} the fundamental transformation, also called, for this reason, the binary 
transformation. 

Another reduction of the fundamental transformation, important from a technical 
point of view, is the Combescure transformation, in which the tangent lines of the 
transformed lattice C(x) are parallel to those of the lattice x. The transformation 
formula reads 



where only the solution Y* of the adjoint linear system (|2.6|) . necessary to build the 
transformation congruence, is needed. 

2.3.3 The Laplace transformations and the geometric meaning of the Hi- 
rota equation 

The Laplace transform Cij(x), i ^ j, of the quadrilateral lattice x is the j-th focal 
lattice of its i-th tangent congruence (see Figure OJ). It is uniquely determined once 
the lattice x is given. The transformation formulas of the lattice points and of the 
r-function read as follows 



d(x)=x-Xi (Y^H^H), 
C*(x)=x-n(X,Y*)(Y*)- 1 H i . 



C(x) = x-Q(X,Y*), 





(2.19) 
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Figure 3: The Laplace transformation Cij 
The superpositions of Laplace transformations satisfy the following identities 

which allow to identify them with the the Schlesinger transformations of the mon- 
odromy theory. 

In the simplest case N = 2 one obtains the so called Laplace sequence of two 
dimensional quadrilateral lattices 

x e = C{ 2 (x), n = C[ 2 (r), = £ 2 i, (.eZ. 

Equations (j2.12J) and (j2.19|) imply that the r-functions of the Laplace sequence satisfy 
the celebrated Hirota equation (the fully discrete Toda system) 

nTiT 2 Ti = (TiT £ )(T 2 n) - (Ti^_i)(T 2 T£ + i). 
2.4 Distinguished integrable reductions 

We will present here basic reductions of the multidimensional quadrilateral lattice. 
The geometric criterion for their integrability is the compatibility with the geometric 
integrability scheme. 

2.4.1 The circular lattices and the Ribaucour congruences 

Quadrilateral lattices Z N — > E M for which each quadrilateral is inscribed in a circle 
are called circular lattices. They are the integrable discrete analogues of submani- 
folds parametrized by curvature coordinates (for example, the orthogonal coordinate 
systems described by the Lame equations (jl.lj) - (J1.2|) ). 

The integrability of circular lattices is the consequence of the fact that, if the three 
"initial" quadrilaterals {x, TiX, TjX, TiTjx}, {x, TiX, T^x, TiT^x}, {x, TjX, T^x, TjT^x}, 
are circular, then also the three new quadrilaterals constructed by adding the vertex 
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Figure 4: The geometric integrability of circular lattices 

circular as well (see Figure 0J). m fact, all the eight vertices belong to 
a sphere, and, in consequence, all the vertices of any K dimensional, K = 2, . . . , N, 
elementary cell belong to a K — 1 dimensional sphere. 

There are various equivalent algebraic descriptions of the circular lattices: 

1. the normalized tangent vectors X{ satisfy the constraint 

Xi ■ TiXj + Xj ■ TjXi = 0, i ^ j; 

2. the scalar function x ■ x : 7j N — ► R satisfies the Laplace equations (|2.1j) of the 
lattice x\ 

3. the functions X° — (x + Tix) ■ Xi : Z N — > R satisfy the same linear system ()2.7|) 
as the normalized tangent vectors X,; 

4. the functions Xi ■ Xi : Z N — > R satisfy equations ()2.9j) and thus can serve as the 
potentials p«. 

The Ribaucour transformation 1Z is the restriction of the fundamental trans- 
formation to the class of circular lattices such that also the "side" quadrilaterals 
{x, TiX, TZ(x), lZ{Tix)} are circular. Again there is no geometric difference between 
the lattice directions and the Ribaucour transformation direction. Moreover, the 
quadrilaterals {x, 1Zi(x), TZ2(x), TZ\ (J^{x)) = TZ% (TZx(x))} are circular as well. In 
consequence, the vertices of the elementary i^-cells, K = 2, . . . , N, of the circular 
lattice and the corresponding vertices of its Ribaucour transform are contained in 
a K dimensional sphere. Finally, for K = N, one obtains a special Z N family of 
iV-dimensional spheres, called the Ribaucour congruence of spheres. 

Algebraically, the Ribaucour transformation needs only a half of the data (nec- 
essary to build the congruence) of the fundamental transformation. The data of the 
vectorial Ribaucour transformation consists of the solution Y* : 7j N — > V*, of the 
linear system (|2.6|) . Then, because of the circularity constraint, Y i: : 7L N — > V given 
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by 

Yi = (Sl(X,Y*) + T i fl(X,Y*)) T X i , 
is a solution of the linear system (|2.7jl . and the constraints 

Q(Y, H) + Sl(X°, Y*) T = 2fl(X, Y*) T x, 
Q(Y, Y*) + Q(Y, Y*) T = 2fl(X, Y*) T fl(X, Y*), 
are admissible. 

We remark that the above constraints have a simple geometric meaning when one 
considers the circular lattices in E M as the stereographic projections of quadrilat- 
eral lattices in the Mobius sphere S M ; i.e., as a special case of quadrilateral lattices 
subjected to quadratic constraints. 

2.4.2 The symmetric lattice 

Given a quadrilateral lattice x with rotation coefficients Qij and potentials given 
by (|2.9jh then the functions Qij, defined by equation 

and called, because of their geometric interpretation, the backward rotation coef- 
ficients, satisfy the Darboux system ()2.8)1 as well. A quadrilateral lattice is called 
symmetric if its forward rotation coefficients Qij are also its backward rotation coef- 
ficients. Again the constraint is compatible with the geometric integrability scheme, 
i.e., it propagates in the construction of the lattice. One can show that a quadrilateral 
lattice is symmetric if and only if its rotation coefficients satisfy the following trilinear 
constraint 

(TiQjijiTjQ^iTkQik) = (TjQi^iTiQt^TkQjk), i,j,k distinct. 

To obtain the corresponding reduction of the fundamental transformation we again 
need only half of the data. Given a solution Y* : I* N — > V*, of the linear system 
(|2.tij) . then, because of the symmetric constraint, Yi : Z, N — > V, defined by 

Y, t = p,,(r 4 r*) T , 

is the solution of the linear system (|2.7|) : notice that, equivalently, we could start 
from Yi. The constraint 

Q(Y,Y*) = H(Y, Y*) T , 

is then admissible and gives a new symmetric lattice. 

There are other multidimensional reductions of the quadrilateral lattice like, for 
example, the D-invariant and Egorov lattices or discrete versions of immersions of 
spaces of constant negative curvature. We remark that the transformations and re- 
ductions discussed above have also a clear interpretation on the level of the analytic 
methods. 
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Figure 5: Elementary quadrilaterals of the isothermic lattice 



3 Integrable discrete surfaces 

In this Section we present some distinguished examples of discrete integrable surfaces. 
Notice that, although the geometric integr ability scheme is meaningless for N = 2, 
however it can be applied indirectly, by considering the discrete surfaces, together 
with their transformations, as sub-lattices of multidimensional lattices. 

We remark also that one can consider integrable evolutions of discrete curves, 
which give equations of the Ablowitz-Ladik hierarchy, and the corresponding inte- 
grable spin chains. 

3.1 Discrete isothermic nets 

An isothermic lattice is a two dimensional circular lattice x : Z 2 — > E M with harmonic 
quadrilaterals,; i.e., given x, T\X and T 2 x, then the point T{T 2 x is the intersection of 
the circle (passing through x, T\X and T 2 x) and the line passing through x and the 
meeting point of the tangents to the circle at T\X and T 2 x (see Figure EJ). Therefore, 
given two discrete curves intersecting in the common vertex Xq, the unique isothermic 
lattice can be found using the above "ruler and compass" construction. 

Algebraically the reduction looks as follows. Any oriented plane in E M can be 
identified with the complex plane C. Given any four complex ponts z±, z 2 , z 3 and z 4 , 
their complex cross-ratio is defined by 



One can show that the cross-ratio is real if and only if the four points are cocircular 
or collinear. In particular, a harmonic quadrilateral with vertices numbered anti- 
clockwise has cross-ratio equal to —1. Therefore, abusing the notation (it can be 
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Figure 6: Asymptotic lattices 



formalized using Clifford algebras), the isothermic lattice is defined by the condition 

q(x,T 1 x,T 1 T 2 x,T 2 x) = -1. 

We remark that the definition of isothermic lattices can be slightly generalized al- 
lowing for the above cross-ratio to be a ratio of two real functions of single discrete 
variables. 

The restriction of the Ribaucour transformation to the class of isothermic lattices, 
named after Darboux who constructed it for isothermic surfaces, has as its data a real 
parameter A and the starting point T>(xq), and can be described as follows. Given the 
elementary quadrilateral {x, T\X, T 2 x, T\T 2 x} of the isothermic lattice, and given the 
point T>(x), then the points T>{T\X) and V(T 2 x) belong to the corresponding planes 
and are constructed from equations 

q(x,V{x),V{T x x),T x x) = A, 
q(x,V(x),V(T 2 x),T 2 x) = -A. 

It turns out that the point T>(T x T 2 x), constructed by the application of the geometric 
integrability scheme, is such that the quadrilateral {V(x), V(T x x), T>(T 2 x), V(T x T 2 x)} 
is harmonic. Moreover, the construction of the Darboux transformation is compatible; 
i.e., the new side quadrilaterals have the correct cross-ratios A and —A. 

There are various integrable reductions of the isothermic lattice, for example the 
constant mean curvature lattice and the minimal lattice. 

3.2 Asymptotic lattices and their reductions 

An asymptotic lattice is a mapping x : Z 2 — > IR 3 such that any point x of the lattice 
is coplanar with its four nearest neighbours T x x, T 2 x, T x ~ x, T 2 ~ x x (see Figure EJ). 
Such a plane is called the tangent plane of the asymptotic lattice in the point x. 

It can be shown that any asymptotic lattice x can be recovered from its suitably 
rescaled normal (to the tangent plane) field N : Z 2 — > M 3 via the discrete analogue 
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of the Lelieuvre formulas 



A 1 x = (T l N)xN, A 2 x = N x (T 2 N). (3.1) 

By the compatibilty of the Lelieuvre formulas, the normal field N satisfies the discrete 
Moutard equation 

T{T 2 N + N = F(TiN + T 2 N), (3.2) 

for some potential F : Z 2 — > M. 

Given a scalar solution 8 of the Moutard equation ()3.2j) . a new solution A4(N) of 
the Moutard equation, with the new potential 

_ mm „ 



m 2 ~ 

can be found via the Moutard transformation equations 

M(TtN) tN = ^-AM{N) T TiN), (3.3) 
lit) 

M(T 2 N) ±N = ^-AM{N) ± T 2 N). (3.4) 

Now, via the Lelieuvre formulas (|3.1jl . one can construct a new asymptotic lattice 
M.(x) = x±Ai (N) x N. The lines connecting corresponding points of the asymptotic 
lattices x and A4(x) are tangent to both lattices. Such a Z 2 -family of lines in IR 3 is 
called Weingarten (or W for short) congruence. Notice that this is not a congruence 
as considered in Section 12. 3. 11 

Various integrable reductions of asymptotic lattices are known in the literature: 
pseudospherical lattices, asymptotic Bianchi lattices and isothermally-asymptotic (or 
Fubini-Ragazzi) lattices, discrete (proper and improper) afline spheres. 

Formally, the Moutard transformation is a reduction of the (projective version of 
the) fundamental transformation for the Moutard reduction of the Laplace equation. 
However the geometric relation between asymptotic lattices and quadrilateral lattices 
is more subtle and the geometric scenery of this connection is the line geometry of 
Pliicker. Straight lines in M. 3 C P 3 are considered there as points of the so called 
Pliicker quadric Qp C P 5 . A discrete asymptotic net in P 3 , viewed as the envelope 
of its tangent planes, corresponds to a congruence of isotropic lines in Qp, whose 
focal lattices represent the asymptotic directions. The discrete W- congruences are 
represented by two dimensional quadrilateral lattices in the Pliicker quadric. 

3.3 The Koenigs lattice 

A two dimensional quadrilateral lattice x : Z 2 — ► P is called a Koenigs lattice 
if, for every point x of the lattice, the six points x±i, TiX±%, T 2 a? ±1 , % = 1,2 of its 
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Figure 7: The Koenigs lattice 
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Laplace transforms (see Section l2.3.3j) belong to a conic (see Figure[7|). The nonlinear 
constraint in definition of the Koenigs lattice can be linearized, with the help of the 
Pascal "mystic hexagon" theorem, to the form that the line passing through x and 
T1T2X, the line passing through x\ and T^X-i, and the line passing through X-\ and 
T 2 2 x\ intersect in a point. 

Algebraically, the geometric Koenigs lattice condition means that the Laplace 
equation of the lattice in homogeneous coordinates x : Z 2 — > can be gauged 

into the form 



It turns out that, if N is a solution of the Moutard equation (|3.2|) . then x = 
T\N + T2N satisfies the Koenigs lattice equation. Therefore, the algebraic theory of 
the discrete Koenigs lattice equation (|3.5j) . its (Koenigs) transformation and the per- 
mutability of the superpositions of such transformations is based on the corresponding 
theory for the Moutard equation ([3.2)1 . 

Geometrically, the Koenigs lattices are selected from the quadrilateral lattices 
as follows. Given a two dimensional quadrilateral lattice x : Z 2 — > ¥ M and given 
a congruence I with lines passing through the corresponding points of the lattice. 
Denote by y i = Tf l l Pi I, i — 1, 2, points of the focal lattices of the congruence. For 
every line I, denote by 1 the unique projective involution exchanging y i with Tiy^ If, 
for every congruence I, the lattice K(x) : Z 2 — > P M , with points JC(x) = i(x), is a 
quadrilateral lattice, then the lattice a? is a Koenigs lattice. The above construction 
gives also the corresponding reduction of the fundamental transformation. 

A distinguished reduction of the Koenigs lattice is the quadrilateral Bianchi lat- 
tice. The natural continuous limit of the corresponding equation is equivalent to the 
Bianchi (or hyperbolic Ernst) system describing the interaction of planar gravitational 
waves. 

3.4 Discrete two dimensional Schrodinger equation 

In the previous secions we have discussed examples of integrable discrete geometries 
described by equations of hyperbolic type. Below we present some results associated 
with the elliptic case; it is remarkable that the QL provides a way to connect these 
two subjects. 

Consider a solution N : Z 2 — > IR 3 of the general self-adjoint 5-point scheme on 
the star of the Z 2 lattice 



T^x + x = Ti(Fx) + T 2 (Fx). 



(3.5) 



aTiiV + T 1 -1 (aiV) + bT 2 N + T 2 \bN) — cN = 



(3.6) 



then the lattice x : Z 2 



M. 3 obtained by the Lelieuvre type formulas 



A ± x 



(T^b)N x T^N, 



A 2 x 



(T^a)N x T^N, 



(3.7) 
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is a quadrilateral lattice having N as normal (to the planes of elementary quadrilat- 
erals) vector field. 

The following gauge-equivalent form of equation (j3.6j) 

1 T 1 ij + T 1 -\^-^) + 7 ^-T 2 ij + T 2 -\^-ij)-q^ = 0, (3.8) 



T X F r 1 v Tir r/ T 2 r r z v T 2 r 

an integrable discretization of the Schrodinger equation 

d 2 ip d 2 ip 
oxf OX2 

is also the Lax operator associated with an integrable generalization of the Toda law 
to the square lattice. 

The 5-point scheme (jH.fi j) is also a distinguished illustrative example of the sub- 
lattice theory. Indeed it can be obtained restricting to the even sub-lattice 7? e the 
discrete Cauchy-Riemann equations 

T 1 T 2 <t>-(t> = iG{T l( l>-T 2 <l>). (3.9) 

Because of the equivalence (on the discrete level!) between equation (j3.9j) and the 
discrete Moutard equation (j3.2j) . the 5-point scheme (j3.6|) inherits integrability prop- 
erties (Darboux-type transformations, superposition formulas, analytic methods of 
solution) from the corresponding (and simpler) integrability properties of the discrete 
Moutard equation. 
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